Theory of non-equilibrium electronic Mach-Zehnder interferometer 
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We develop a theoretical description of interaction-induced phenomena in an electronic Mach- 
Zehnder interferometer formed by integer quantum Hall edge states (with v = 1 and 2 channels) out 
of equilibrium. Using the non-equilibrium functional bosonization framework, we derive an effective 
action which contains all the physics of the problem. We apply the theory to the model of a short- 
range interaction and to a more realistic case of long-range Coulomb interaction. The theory takes 
into account interaction-induced effects of dispersion of plasmons, charging, and decoherence. In 
the case of long-range interaction we find a good agreement between our theoretical results for the 
visibility of Aharonov-Bohm oscillations and experimental data. 
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I. INTRODUCTION 

Many recent experiments studied transport through an 
electronic analog of Mach-Zehnder interferometer (MZI) 
built on edge states in the quantum Hall regime 1 - - — . 
These experiments show strong Aharonov-Bohm oscilla- 
tions, which is a manifestation of quantum interference 
of electrons propagating along the arms of the interfer- 
ometer. One of remarkable experimental observation is a 
lobe- type structure in the dependence of visibility on bias 
voltage. More precisely, the visibility strongly depends 
on voltage, showing decaying oscillations ( "lobes" ) char- 
acterized by certain energy scale. Such structure can not 
be explained within a model of non-interacting electrons 
(which would predict a constant visibility) and thus re- 
sults from the electron-electron interaction. Therefore, 
the experiments on Mach-Zehnder interferometers ex- 
hibit the physics resulting from an interplay of a quantum 
interference and the Coulomb interaction under strongly 
non-equilibrium conditions. Development of a theory of 
such phenomena is a challenging task. 

Various aspects of the relevant physics have been ad- 
dressed in earlier works. Dcphasing of quantum inter- 
ference in Aharonov-Bohm rings, interferometers and re- 
lated phenomena have been studied in many works for 
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and non-equilibrium^— setups. The 
importance of electron-electron interaction for the emer- 
gence of the lobe structure has been emphasized in 
Refs. iHl In Ref. [H the influence of long-range (1/r) 
character of the Coulomb interaction (leading to a dis- 
persion of the plasmon mode) was analyzed. In the works 
12514271 a simplified model with a quantum-dot-type treat- 
ment of the interaction was solved. 

The goal of the present article is to present a systematic 
theory of transport in a non-equilibrium quantum Hall 
Mach-Zehnder interferometer and to confront its predic- 



tions with the experiment. This theory formulated in 
terms of functional-bosonization Keldysh action in Sec.lHl 
takes into account all effects of the electron-electron in- 
teraction, including formation and characteristics of the 
lobe structure, as well as the dephasing. In Sec. IIIII we 
apply the theory to the model of short-range interaction. 
In this way we reproduce the earlier results on the lobe 
structure^ and also find the suppression of the interfer- 
ence signal due to dcphasing. While showing some sim- 
ilarity to experimental observations, the results of the 
short-range-interaction model contradict to the experi- 
ment in several crucial aspects. This motivates us to 
explore the more realistic case of 1/r Coulomb interac- 
tion in Sec. IIV1 Using our general formalism, we analyze 
the cases of v = 1 and v = 2 edge modes. The obtained 
parametric dependences of the interference signal are in 
good agreement with experiments, although in the case 
of v = 1 mode some discrepancies remain. We summa- 
rize our findings and discuss further research directions 
in Sec. IVl 



II. GENERAL FRAMEWORK 



Model and the functional bosonization 



We consider a theoretical model of the Mach-Zehnder 
interferometer, realized with edge states in the quan- 
tum Hall regime at filling fraction v = 1 and v = 2, as 
schematically shown in Figs. [T] and The outer channels 
propagating along different arms (index: ±) of the inter- 
ferometer are coupled by means of two quantum point 
contacts (QPCs), located at points xf and xf. Each 
QPC can be generally described by the unitary 2x2 
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scattering matrix 



1,2 



(1) 



where the transmission and reflection coeflicients, Tj and 
Tj are assumed to be real. In case of v = 2 this scatter- 
ing matrix relates incoming modes of two outer channels 
(in up/down arms) with the outgoing ones. Electrons in 
the inner channels are assumed to propagate through the 
QPCs without scattering. 

The arms of the interferometer can be different, 
thereby generally ^ xj . The coordinates and re- 
fer to the points where the MZI is connected to the drain 
and source reservoirs. Specifying the non-equilibrium 
boundary conditions, we assume that only one single 
channel is biased at the chemical potential eV, while all 
others are grounded. In case of v = 2 it will be the outer 
upper channel, as shown in Fig. [I] In the experiments^^ 
such situation is realized with the use of an extra QPC 
which splits the incoming inner and outer channels [sec 
Fig. U], so that these channels originate from different 
reservoirs. 

To set the stage, we start by considering the simplest 
situation, when Tj = 1 and Tj = 0, so that outer chan- 
nels in the two arms of the interferometer are completely 
decoupled from each other. In this case we model the 
system by a set of interacting chiral fermions with the 
action S = Sq + S- m t , where 

S a = -iv^2 / dtdxip^(x)d x ip a (x), (2) 

a J 

Sint = 7}^2 J dt dxdx' p a (x)U (x - x')p a (x') . 




QPC2 



FIG. 1: Scheme of an electronic Mach-Zehnder interferometer 
built on quantum-Hall edge states at filling factor v = 2. 
Quantum point contact (QPC1 and QPC2) characterized by 
transmission amplitudes Tin) are used to partially mix the 
outer edge channels. All Ohmic contacts are grounded, except 
for the source terminal 5*3 which is kept at voltage V. The 
current is measured in the drain terminal D%. The QPCO is 
pinched in such a way that the inner channel is completely 
reflected while the outer one is fully transmitted. 
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FIG. 2: Top: Schematics of a Mach-Zehnder interferometer. 
Two quantum point contacts are characterized by scattering 
matrices Si and S2. The electron-electron interaction is de- 
coupled via a Hubbard-Stratonovich fields <p + and ip~ on two 
arms of the interferometer. Bottom: Coordinates on two arms 
of the interferometer. 

Here ip a is a Grassmann field of the chiral fermion in the 
arm a = ±. In case of v = 2 it has the vector structure 
i'a = {i } ioni } 2a) T : the subscript p = 1,2 denoting the 
outer and the inner channels, respectively. The potential 
Uq describes the bare interaction within the interferom- 
eter arm between the chiral electron densities 

p a {x) = ipl(x + 0)tp a (x), 

and v denotes the drift velocity in the edge. 

In what follows we are going to use the Keldysh ver- 
sion of the functional bosonization^. Let us recall the 
basic ideas of this construction. We decouple the inter- 
action term 5i n t using the Hubbard-Stratonovich trans- 
formation with a field <p and double the number of fields; 
ip = (ipf,ipb) T (and the same for <p), where tpftb) denote 
the Grassmann fields residing at the forward (backward) 
branch of the Keldysh contour. These steps lead us to 
the action 



A 



E 

a 

1 

2 



dtdxipa(id t + ivd x - ip a )ip 



dtdxdx' (p a U 1 (x — x')ip 



Ck 



(3) 



where we assume the implicit summation over the chan- 
nel indices. 

The minimal coupling between fermionic and bosonic 
degrees of freedom in ID geometry can be eliminated by a 
local gauge transformation, ip a — > ^ a e ieo , if one requires 
that 



(d t +vd x )G c 



-fa- 



(4) 



One has to resolve this differential equation by taking the 
proper structure of the Keldysh theory into account, 



e 6 



(5) 
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Here we have denoted £ = (x, t) and omitted the 
arm/channel indicies for brevity. At zero temperature 
elements of the bare particle-hole propagator Dq in the 
space-time representation are given by the relations 



D, 



T/T 



(0 



v 1 n^(t ~ x/v), 

9(Tt)D<(0 + 9(±t)D>(0, (6) 



where n^(i) = —i/2Tr(t =F ia/v), 9(t) is the Heaviside 
theta-function and a is a short (ultraviolet) cutoff scale 
that is of the order of magnetic length lg. After this 
gauge transformation the Green function of interacting 
electron acquires the form 

G > (a,6)=ff > (6-6)(e^ (6) e-^fe^ . (7) 

Here g > (£) is the free zero-temperature Green function, 



9 > (0 



-, — i(t — xfv)eV 



27TV (x/v — t 4- ia/v) 



(8) 



that should be understood as a diagonal 4x4 (two arms, 
two channels) matrix; for the sake of generality, we as- 
sume that at v = 2 channels are biased by distinct volt- 
ages Vp , where p = 1, 2 is the channel index, and a = ± 
is the arm index. The average over phases 0(£) in Eq. ([7]) 
is performed with the Gaussian action 

Mf) = \^^ x u z (p- ^n^-^azpo, (9) 

which is a special property of the ID geometry (Larkin- 
Dzyaloshinskii theorem). Written in the symbolic form, 
this expression implies the summation over the Kcldysh 
indices (determining the vector structure of 0) and con- 
volution in the space/time domain, see Eq. ([5]). The ma- 
trix n is the polarization operator (a, a' denoting the 
Kcldysh indices) 

n aa '(£) = -ig aa '(0g a ' a (-0 - (2™)-^)^', (10) 

with the last term accounting for the static compress- 
ibility of edge channels. The linear in <p term in the 
action © describes the response of the system to the ex- 
ternal non-equilibrium charge p l = (po, po) injected from 
the leads 



Po (Z) = (2*v)- 1 eV(t-x/v). 



(11) 



It is convenient to perform the Keldysh rotation (see e.g., 
Refl30l) , introducing the "classical" and "quantum" com- 
ponents of the fields 



¥>c{q) = iff ± <Pb)/2, 



(12) 



Then the polarization operator in the frequency- 
momentum representation acquires the Keldysh struc- 
ture with 



n 



R,A 



<1 



2tt qv — to =F iO 
n A - = (U R -U A )coth(f3cj/2). 



(13) 



Minimizing the action (|9|), 8Ao{(f) / S(p q = 0, one obtains 
the linear equations for the mean value of a classical elec- 
trostatic field (p c and an induced charge pi at the edge, 



u o (<?)</>c = Pi + po, 
pi = n_R(g, uj)ip c . 



(14) 



Below we concentrate on the experimental situation with 
a dc applied voltage which introduces a homogeneous 
charge, so that we can take the limit uj = 0, q ~ l/£, 
with L being the system size. 

Evaluation of the Green function Q is now reduced to 
a Gaussian integral: 

G>(£,0 = 9> tf-O exp {i[0(O - 0(O] -J>(Z- £')} > 

(15) 

where the average phase is induced by the mean clas- 
sical electrostatic field 



e(0 



dx' tp c {x 1 , t + (x' — x) /v) , 



(16) 



and J is the correlation function of the Gaussian phase 
fluctuations (50 around the above mean value: 



1 



J aa (0 = 7;([se a (0-se a (o)Y)ao 



(17) 



The phase-phase correlation function can be expressed 
in terms of the bare particle-hole propagator D and the 
effective interaction 



{w t ) = -u = -{a z Uv l -ny 1 . 



(18) 



Using the linear relation ([5]) between the phase and the 
electrostatic potential, one obtains 



where 



J aa (£) = V aa (£) - V aa (0), 



(19) 



(20) 



The correlation function J aa (£) can be now easily eval- 
uated with the use of (ui, q) representation. The result 
reads 



(21) 



J aa (0 = v Jp (0 - J a F a (0 



where we decomposed J(£) into plasmon (P) and free (F) 
parts, 



J P/F 



dq 

with the plasmon velocity 

u P (q) = v + vU (q) /{2tt) 



(22) 



(23) 



and the drift velocity up = v, respectively. Equation 
(f2Tj) is the zero temperature result for the "greater" func- 
tion. The "lesser" correlation functions satisfies J < (£) = 
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[J 2 * (£)]*, and the (anti-) time-ordered correlators can be 
reconstructed with the use of basic definitions, akin to 
Eq. ©. 



Equation (|2"2"j) yields for Jp\ 



In 



x — vt + ia 



i a 



(24) 



Thus, in the case v = 1 the bare electron pole in Eq. is 
canceled by the free contribution to J(£) and the electron 
motion is determined solely by the plasmon contribution, 
whereas for v — 2 both plasmon and free (neutral) modes 
do influence the behavior of electronic propagator. The 
plasmon contribution Jp(£) depends on the specific form 
of the electron-electron interaction through the momen- 
tum dispersion in the velocity up(q). We consider two 
models of short-range, Uq(x) oc 5(x), and long-range, 
Uq(x) oc 1/\x\, interactions in the following Sections ITTA 
and El 



Keldysh action of the Mach-Zehnder 
interferometer 



In this section wc generalize the discussion of the 
functional bosonization to include electron scattering 
at QPCs which are crucial elements of the MZI lay- 
out. In this way we construct the most general Keldysh 
action A{<p) of the MZI, which is applicable for arbi- 
trary scattering matrices characterizing two QPCs. The 
Keldysh action can be expressed in terms of a total single- 
particle time-dependent interferometer scattering matrix 
S(t,t', [ip]). It describes scattering of electrons at QPCs 
as well as their propagation along the arms of the inter- 
ferometer in the fluctuating Hubbard-Stratonovich field 
(p. The 5-matrix is thus non-local in time and takes 
different values (Sf = S[tpf] and Sb = S[ipb) ) on the for- 
ward (backward) branch of the Keldysh contour, since 
the fields ipb(x,t) and ipf(x,t) generally differ. The ac- 
tion A(ip) has the form 



iA(ip,x) 



lndet 



I + (Sle^S f - 1)/ 



— iip n (p 



(25) 



Here <p — (tpcPq) 1 is assumed to be written in the 
Keldysh basis. We have also introduced the elec- 
tron distribution function in the source reservoirs, / = 
diag(/ + , /~) (in the case v = 2 each f ± is itself a diago- 
nal matrix in the channel space with the components fp), 
and auxiliary "counting fields" \ = diag(x + ,x~) in the 
drains which enable us to find the number of transferred 
electrons. The determinant in Eq. (|25p is evaluated with 
respect to time, arm, and channel indices. Since in this 
paper we restrict ourselves to the case of zero tempera- 
ture, the distribution functions in time domain read 



ft(t,t') 



' 2m(t - t' + iO) ' 



(26) 



where correspond to voltages applied to the different 
channels. The polarization operator II has the form 



If 



o n. 4 
o 



(27) 



We now specify the total scattering matrix S for the 
MZI. It can be decomposed into the scattering matrices 
of the quantum point contacts (QPC), S\ and S2, and the 
transfer matrices accounting for the propagation from 
("source") to x\ (Tio), from x\ to x-i (T21) and from X2 
to L (T L2 ) ("drain") (see Fig. EJ, 



S — TL2S2T21S1T' 



10 ■ 



(28) 



with the scattering matrices of the contacts defined by 
Eq. |T]). Since the Hubbard-Stratonovich field is diagonal 
in the channel space, the transfer matrix can be written 
as 

r i4 (t,f) = Ajiit^eW), (29) 
where the shift operator 

Af i (t,t')=d(t-t'-(xf-xf)/v) (30) 
describes the free propagation, and 

6%(t) = -v- 1 l dx'^(x\t+(x'-xf)/v) (31) 



is the electron phase induced by the Hubbard- 
Stratonovich field during the propagation of an electron 
from Xi to Xj. The weak variation of magnetic field 
gives rise to an additional Aharonov-Bohm phase $ = 
$ + — $_, which we incorporate by Q^i ~^ ^21 + ^ >± - Def- 
inition (|28| gives both 5/ and Sb depending on weather 
iff or ifb is used to reconstruct the electron phase in 
Eq. (|3"Tj) . Wc note, that the phase df/ b > considered as 
the function of (xj,t), satisfies Eq. ((4]). However, this 
phase is different from 6y/ h , because the relation (|3Tj) 
does not contain information about Keldysh distribution 
functions, as opposed to Eq. (0). 

The cumulant generating function of the MZI can be 
now expressed as a functional integral over ip, 



Z(x)= / £V// 6 0M) exp{iA(ip,x)}- 



(32) 



In particular, it gives the mean number of transferred 
charges to the upper/lower drains, 



-id x ± lnZ(x)lx=o- 



(33) 



A detailed derivation of the action (|25|) will be pub- 
lished elsewhere^. This action bears connection with 
the solution of the problem of full counting statistic s) 31 ' 32 . 
Related Keldysh actions of the determinant structure 
were found for the problems of a local scatterer in Ref . 
and of a quantum wire in Ref. l2ll . One can also show that 






FIG. 3: Example of paths contributing to the matrix element R — Q++. The paths (a) and (b) give direct contributions, while 
the path (c) encloses the magnetic flux $ and describes the interference process. The corresponding analytical expressions read: 



TMrrnr«r{ ; (b) 
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(c) rirariraryr'e^e 



in the limit of zero tunneling, n = 0, the full action (f25j) 
is reduced to the Gaussian one, given by Eq. ©. In this 
case the linear in ip term and the classical Gaussian fluc- 
tuations of <p, described by II k [sec Eg. ([T3")) ]. happen to 
be the first two terms in the expansion of the functional 
determinant of the action (j2"5"]) . while the higher order 
terms in tp vanish. 

In general, the integral over the Hubbard-Stratonovich 
field with the determinant action A(p) can not be eval- 
uated analytically. In the following, we proceed with the 
analysis of the MZI in the weak tunneling limit, tj <C 1. 
Since in the absence of tunneling the MZI is described by 
the Gaussian action (|9]), we can introduce the tunneling 
action At{<p), so that A(p) = Ao{p) + At{p), where the 
expansion of At in terms of Tj starts from the terms of or- 
der of 0(t 2 ). In Appendix A we show that the tunneling 
action At can be obtained by a proper regularization of 
the functional determinant appearing in Eq. (|25p . yield- 
ing 

1 + (e l ^Sl(x)S f ( X )e-^f - l) f] . (34) 



iAt = lndct 
Here 

S s { X ) = ^ /2 S s e-^l\ sl( X ) = e-^ 2 Ste^, (35) 



and 



dx 



ip fb (x',t + x'/v) 



(36) 



is the phase collected along the way from to L^ 1 with- 
out tunneling. In other words, ■0 ± is equal to Of^ with 



= and 



- T± 



L . New distribution functions are 



obtained from the original distribution functions in the 
sources (|26[) with the use of a gauge transformation, 



f(t,t') = e^fit^e-^ 
- f>$f, a + f<$t, a > 



(37) 
(38) 



where we defined the projectors on the hole and particle 
subspaces 



f>(t-t') =±l/2m(t-t' TiO) 



(39) 



(see also Appendix A). According to our general nota- 
tional conventions, Eq. (|3"5f has to be understood as a 
convolution in time domain. 



To proceed with the evaluation of the tunneling action, 
we introduce a matrix (in the time, arm, and channel 
space) 



Q = e^sl(x)S f (x)< 



-iijj S _ 



R T 
T R' 



(40) 



The second equation in Eq. (|4TJ| introduces a block de- 
composition of Q in the arm space. In the absence of 
tunneling Q(t,t') = S(t — t') and At vanishes. Since the 
tunneling affects only outer channels, it is sufficient to 
keep only the corresponding elements of the matrix Q. 
Therefore, from now on Q is the matrix in the arm and 
time space only, both for v = 1 and v = 2 setups. 

To find Q in the presence of tunneling we define the 
"hopping" operators at each QPC, 



(t,t') 



,f/b 



m-t'-{ X +- X -)/v)e-<l f,b ^'\ 

(41) 

where 0"'/^ b (t) is the phase picked up by the electron 
along the path from Xj to the drain a, as given by the 
relation (|3"Tj) . We can further formulate a set of simple 
rules that enable one to express the matrix elements Q a p 
in terms of these operators: 

(i) Each matrix element Q a /3 is associated with a num- 
ber of possible continuous paths that originate from 
the source terminal a, go to any of the drain termi- 
nals and return back to the source /3 (see Fig. [3]) . 
The forward and backward paths are allowed to go 
through the different interferometer arms. In the 
latter case the associated contribution corresponds 
to the interference process. Each element Q a p is 
thus a sum of eight contributions. 

(ii) Each transmission via QPC j gives the factor Tj. 
Each reflection gives the factor ±irj along the for- 
ward/backward path respectively. 

(iii) Each tunneling process from arm + to — which 
happens at the QPC j is associated with the am- 
plitude t{' b e T%x ^ 2 , where x — X+ — X-- The in- 
verse tunneling process is described by the Hermi- 
tian conjugate of this operator. The product of 
such "hopping" operators must be path-ordered. 

(iv) contributions associated with the interference pro- 
cess carry the extra Aharonov-Bohm phase e ±z4> . 
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As an example, three (out of eight) contributions to the 
matrix element R = Q++ are shown in Fig. [3] The full 
analytical expressions for Q a p is presented in Appendix 
B. The next step is to use the expansion 

iA t = Trf(Q-l)/-i((Q-l)/) 2 ) +G(r 3 ). (42) 

Some important technical details of the evaluation of At 
with the use of this formula are given in Appendix C. We 
show there that the second-order approximation in Tj for 
the tunneling action At has a form 



iA t 



= J2 T ^ 

ij 

x / di li2 e ie ^" tl )n ij -(t 1 ,t 2 )e- ie( ^' t2) . 



(43) 



Here the summation is performed with respect to QPC 
indices (i,j = 1,2) and the time integration goes along 
the Keldysh contour. The generalized polarization oper- 
ators for the direct terms read 



n 



T/T 



(ti.ta) = - />,<(*! - *a)/<,>(*2 - *i) , 



(44) 



and the polarization operators of the interference terms 
are given by 



n 



T/T 
12 



{tl,t 2 ) 



n&(ti,t 2 ) 

I^/ f (*l,*2) 

n|(ti,t 2 ) 



-e- a ft,>(ti-t2-{x+-xt)/v) 

x/>,< (t 2 - *i - (»2 - Xi)/v) , 
nf/ T (t lj t a )e= H x ) (45) 
- e »/^<(t 1 -t 2 -(4-^)/«) 

x/<,> (*2 - *1 - - ^2 )/ u ) . 

n 2 n / f (t 1 ,t 2 )e ± ^. (46) 



They are built up from the electron and hole distribution 
functions 



ft(t-t') = e- ieV± ^Mt-t'). 



(47) 



The polarization operators are dressed by relative gauge 
factors Q(xi,t) = Q + (x^,t) — Q~(x~,t), where the 
phases at each arms ± arc linked to the corresponding 
HS-field ( / 3 ± by relation ([5]), thereby taking the interac- 
tion into account. 

The structure of the polarization operators can be in- 
terpreted as follows. The component II < oc /</> cor- 
responds to a particle excitation in the arm + and a 
hole excitation in the arm — . This implies the charge 
transfer from + to — , which explains the "counting" fac- 
tor e~ %x . Similar considerations explain the structure of 
n > . Finally, the components II T / T at \ — are related 
to IK via the conventional Kcldysh-formalism relations, 



C. Gaussian approximation 

The Keldysh action enables us to calculate the tunnel- 
ing current of the interferometer. Specifically, let us con- 
sider the current of particles transported from the upper 
(+) to the lower interferometer arm (— ). This current 
can be obtained as derivative of the cumulant generating 
function Z(x) with respect to \- ■ We express the current 
as difference of up- and down rates, I = ^2^(1^ — iy), 
with 



to J-oo 



. e te i '- b (x i ,t 1 ) e -ie b ^(x j ,t 2 ) 



A +A t (x=0) 



(48) 



where we introduced the measurement time to . The aver- 
age over phases is to be done with the full action. How- 
ever, to the leading approximation, we can neglect the 
tunneling action and perform the phase average with the 
Gaussian action Aq. The corresponding correlation func- 
tions of phase variables are given by Eq. (|21[) . Assuming 
that no interaction is present between the different inter- 
ferometer arms, wc find 



(49) 



dtG+'*(xt -x+ t)G~&(xr -xr,-t), 



where the full Green functions are given by Eq. (|15p . 
Equation (|49|) was the starting point of the discussion 
in the works by Chalker, Gefen and Veillette^ and Lev- 
kivskyi and Sukhorukov^. We go beyond it by including 
into consideration the (non-Gaussian) tunneling action 
that controls the non-equilibrium dephasing. 



D. Non-Gaussian effects 

Let us now improve the result for the tunneling cur- 
rent (|49p by taking into account the non-Gaussian na- 
ture of phase fluctuations, which arc described by the 
tunneling part of the action At- These fluctuations are 
responsible for the intrinsic dephasing in the MZI by the 
non-equilibrium shot noise generated in the QPCs. To 
be specific, we evaluate which requires calculation of 
the correlation function 

^ e iQ f (x l ,t 1 ) e -iB h (x J ,t 2 )\^ _ j -p^iA^+iAo+iAt^o^ 

(50) 

where we have introduced a linear-in-</? term in the ac- 
tion, 

Af^tuh,^}) = e f (x i ,t 1 )-e b {x j ,t 2 ) 

= J2j dtdx(Q a Yj a . (51) 
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Here — {J^'* ,S are source terms acting on 
both branches of the Keldysh contour, 

J ±J {x,t) = ±6{x-xf)S(t-t 1 ), 

J^ b (x,t) = ^6(x-xf)6(t-t 2 ). (52) 

In order to account for the non-Gaussian effects, we 
proceed with a real-time instanton approach developed 
in Ref. [H for the problem of tunneling spectroscopy of 
a biased quantum wire with a scatterer. We look for the 
optimal (saddle-point) trajectory ip"(x,t), which mini- 
mizes the total action Atot — Aq + At + Afj at % = 0. In 
the weak tunneling limit, Tj <C 1, one can find the instan- 
ton trajectory ip* approximately by minimizing only the 
quadratic part of the action, i.e. Aq + Afj. This yields 
for the stationary point 

6 Q [^] = -VJ a , (53) 

where the particle-hole propagator is defined in Eq. (f20| . 
One can show that taking into account corrections of 
order r 2 which follow from the exact non-linear equation 
of motion would lead to a contribution of order r 4 to the 
tunneling action. Thus, such corrections are negligible 
within the accuracy r 2 of our calculation. 

For simplicity, we do not take into account at this stage 
charging effects on the interferometer arms. These effects 
will be restored at the end of the calculation as constant- 
in-time phase shifts evaluated according to Eq. ([To]) . 

Following this procedure, we find for the tunneling ac- 
tion in the saddle-point approximation 

'A:;,, = £7*7, [ dt[ t2 e i& l^m u (t[,t' 2 )e- ie n^. 

(54) 

with the phases 

ierw = j fa (xt-x+, ti-t)-J ba {x}-xt, t 2 -t) 

+ J fa (xr -x-^-t)- J ba {xj -x^,t 2 -t) 

(55) 




FIG. 4: Diagrammatic series for the tunneling current corre- 
sponding to the real-time instanton configuration ()55|) . Black 
wavy lines, representing the exponent of the phase correlation 
function, exp{ — J(£)}, are taken into account by the saddle- 
point approximation, while the red ones originate from the 
quantum fluctuations around the instanton. 



(a is the Keldysh index). In the stationary situation this 
action should become a function of time difference only, 
A* t ^{t\ —t 2 ). We have also taken into account quantum 
fluctuations around the instanton trajectory ip". One 
can prove that, in analogy with Ref. |22| . they result in 
renormalization of the bare polarization operator Tl^i to 
the dressed one given by 

fLit'(t u t 2 ) = K?( tl ,t 2 ) 

x exp{-J aa '(t 1 -t 2 ,x+ - x-)} 

x exp{-J a ' a (i 2 - t±,xl - x^" )}. (56) 

The tunneling current can be now evaluated by sub- 
stituting the stationary phase O* into relations (|4"5)l and 
(|50p . The diagrammatic representation resulting from 
expansion of the exponential of the tunneling action 
evaluated in the saddle-point approximation is shown in 
Fig. 2J If we retain only the first term of this expansion 
(i.e. neglect the tunneling action), we reproduce the re- 
sult (j49|) of the previous section. The non-Gaussian shot 
noise modifies this result by the additional factor e lAt 
under the time integral, 

/| = er^v 2 J dt'exp{i^(t')} (57) 
x G + ^{xf - xf,t')G->$(xJ - xr, -t'). 

The full Green's function G + here contains the extra 
phase shift, which takes into account the non-equilibrium 
charge effects. It has to be found in accordance with 
Eqs. (fl4"|) and (|I6j) . Solving mean field equations (|I4[) in 
the dc limit, i.e. at to = 0, q ~ l/£, with the use of 
po = eV/(2irv) as the injected non-equilibrium charge, 
one gets the mean electrostatic potential on the upper 
arm of the interferometer, 

fe+ = ""( 1 -MT7I)) el ' l58) 

and the non-equilibrium phase shift 

et(x) = u- 1 ( 1 -l)eVx. (59) 
\up(l/L) vj 

Using the expression (|57p , we can now derive a general 
result for the differential conductance Q{V) = dl/dV. 
Specifically, the conductance is obtained as a sum of two 
contributions: Gq(V) which is independent on the en- 
closed flux $ and another one, G<s>(V), which is sensitive 
to 

2 

G{V) = ^— (r 2 Ii + t 2 1 2 + 2nr 2 Rc (Z$e 4 *)) , (60) 
where 

l k ?t Attv 2 J dtte- teVt Im(G{;(0,t)) 2 cxp{-lmAl kk (t)} , 

1$ ~ Anv 2 J dt(t- At) e- ieV{t - At) exp{-Im A* t 12 (t)} 

x Im(G>(l+,t)G<(-r,-t)). (61) 
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In this expression I 



are the lengths of the 



upper and lower arms of the MZI, and 



At = 



d 



d{eV) 



e?04)-e c + (zr) +i+/v 



(62) 



is the delay time related to the charging effects discussed 
in Sec. Ill Al The explicit form of the delay time At de- 
pends on specific model of the interaction. When deriv- 
ing the result (|5T|) . we have neglected the real part of 
the tunneling action, Re (it leads to small corrections 
only that do not affect the result in any essential way). 



III. SHORT-RANGE INTERACTION 

In this section we analyze a model of an interferometer 
with v = 2 channels per interferometer arm and short- 
range (point-like) interaction, U${x — x 1 ) = UqS(x — x'). 
This model was considered previously in Ref. [2J where 
non-equilibrium dephasing effects were discarded. We 
will show that predictions of this model can not be rec- 
onciled with the experiments, which can be traced back 
to the simplified treatment of the interaction. Neverthe- 
less, it is natural to start our analysis from this model 
that serves nicely for illustration of general principles. 

The elementary excitations induced by the short- 
range interaction interaction are dispersiveness collective 
charge and neutral modes with velocities u = vp + Uq/t: 
and v respectively. The correlation functions read 



T > / \ 1 , ( x — ut 

J<(x,t) = - In 

v ' ' 2 \x-vt 



J T ' f (x, t) = 8(t)J* (x, t) + 6(-t)J^(x, t), 



(63) 



(64) 



where t± = t ± iO. It is easy to verify that Eq. (|6"3")l can 
be also written as 



J T ^(x,t) = 6(x)J^(x,t) + 6(-x)J^(x,t) 



(65) 



The last equation is a consequence of the chirality of the 
system. Similar relations hold also for the polarization 
operators Hki , and these relations remain unchanged af- 
ter renormalization. 

Using the relations (|65|) for the correlation functions, 
one can show that the saddle-point phase 0^' a does not 
depend on the Kcldysh index a if k = 2. (We remind 
that, according to our convention, x 2 > x^ .) Analo- 
gously, the renormalizcd polarization operator Iljy does 
not depend on a (a'), if k = 2 (I = 2). Using the above 
relations and taking into account contributions from the 
forward and backward Keldysh contour, we find that the 
only non-vanishing term in the sum (|54p for the tunnel- 
ing action A? t is the one with k = I = 1. In other words, 
the tunneling action in the saddle-point approximation 
has the only contribution oc jf, which means that only 
the noise produced by the first contact contributes to 



dephasing. This remarkable property has a simple phys- 
ical explanation: in chiral ID-system particles possess 
a well-defined prehistory, where the direction in space 
is linked to the direction of time. At zero temperature 
but finite voltage the tunneling of electrons via the first 
QPC is accompanied by a spontaneous emission of non- 
equilibrium plasmons and neutral modes. They transfer 
the shot noise oc Tj 2 to the second QPC, where it affects 
the direct (I22) and interference (/12) current. On the 
other hand, the shot noise generated at the second QPC 
is transferred directly to the drains without affecting the 
particle interference. 

Let us evaluate the tunneling action for the interference 
current 7^, 

*A*,i2 = rf J dtf^eWJfLufo - t' 2 )e-^\ (66) 



(67) 



where = x 2 — xf are the lengths of the interferometer 
arms. Remarkably, this action for the interference cur- 
rent does not depend on the times £1,2, and therefore, the 
current factorizes in the one found in Sec. Ill CI and a de- 
phasing factor. The dominant large-voltage behavior of 
(1551) is determined by the singularities of the polarization 
operator at t[ ~ t' 2 ~ 1/eV, yielding 



iA 



t,12 



dt[P>(e ie ?W- ie * f M-l) 
+P<(e ie * / «- i <(*) -1)], (68) 



where 



P*i= f dttif^t) = ~-9( T eV) 
J u 



\eV\_ 
2tt ' 



(69) 



Thus, for the imaginary part of the action we have 

ImXi2= - 2r 1 2 [P 1 > 1 +P<] (70) 



dt sin 2 



'&f{t)-@\ f {t) 



The only effect of quantum corrections to the saddle- 
point approximation is a renormalization of the tunneling 



strength, r-f 



-2 



tIv/u. 



The final result for the dephasing rate depends on rela- 
tion between the parameters of the problem. In the case 
l + jl~ > u/v, we find 



lmAt >12 = f!-\eV\(l + + l-)(--- 
ir \vu 



(71) 



Therefore, in this regime the dephasing rate is propor- 
tional to the total length of the interferometer. Clearly, 
there is no dephasing when the velocities are equal to 
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FIG. 5: Visibility at v = 2 calculated in the model of strong 
short-range interaction, ti>«, with equal QPCs, f\ = 0.2. 
Blue curve: dephasing is neglected; red curve: dephasing by 
the shot-noise is taken into account. Upper panel: l + /l~ — 
1.15; lower panel: l + /l~ = 1.35. 



each other, which is the case for the non- interacting 
system. If parameters are chosen in such a way that 
l + /l~ < u/v, we get 



i 



\eV\ 



(72) 



Remarkably, we find no dephasing in the case of equal 
arm lengths. We offer the following physical interpreta- 
tion of this result. A tunneling event at the first contact 
is accompanied by emission of plasmons in the upper and 
lower interferometer arms. If the arm lengths are equal, 
the plasmons meet together at the second contact com- 
pletely phase-coherently. In other words, when electrons 
tunnel at the first contact, there will be noise created in 
the upper and lower interferometer arm, and this noise 
is fully correlated. In the case of equal arm lengths it is 
possible to restore the complete phase information at the 
second contact. Thus, dephasing is absent. 

In the case of direct currents, an analogous calculation 
for the long-time limit gives 



ImXiiW 



0. 



(73) 



lmAl 22 (t) = rl-\eVt\. 

7T 



As for the time delay At related to the charging of the up- 
per arm of the interferometer, we obtain, using Eqs. (1141) . 



T6J) and d62]), 



+ /I 1 
At = — - + - 

2 V u v 



(74) 



We now have all the ingredients needed to evaluate 
the visibility of the MZI according to Eqs. ([60]) and (foTj) 
for the differential conductance. Wc define the volt- 
age dependent visibility as the ratio of the amplitude of 
Aharonov-Bohm oscillations in conductance to its mean 
value: 

V G (V) = (max(£ $ 0O) -min(0*(y)))/2G o (tO- (75) 

We observe that the direct conductance £/o(^0 does not 
depend on voltage despite a non-zero action A* t22 : the 
latter affects each of the tunneling rates and 1^% but 
does not change their difference. Therefore, the visibility 
can be expressed via the integral 2$ given by the rela- 
tion (|6"Tj) as 



V G (V) = 



2nr 2 



(76) 



In what follows we consider the limit of strong interac- 
tion, «»!). In this case we represent 1$ in the form 



1$ = 2^ exp{- 
where 1^ is simplified to 



r(0) 



f dt, 
J -(i-At)e- 



ieV(t-At) 



(77) 



(78) 



x Im 



(t + i0)(t - ^ + i0)i(t - £p + iO) 



Two examples of the visibility plots in the strong- 
interaction limit and equal transmissions at both QPCs 
are shown in Fig. [5j where we have denoted / = (l + + 
l~)/2. 

While the obtained voltage dependences of visibility 
does show an oscillatory structure similar to experimen- 
tally observed lobes, there are strong differences. Most 
importantly, in the limit of equal arms, which is predom- 
inantly the experimental situation, the visibility calcu- 
lated in the framework of short-range interaction model 
oscillates as Vg oc cos(eVl/v) without any decay, since 
the non-equilibrium dephasing rate vanishes according to 
the relation (|T2"|) whatever strong the interaction is. The 
situation gets even worse when we try to apply the same 
model to describe experimental data on MZI at v — 1. 
Specifically, the model predicts then a constant visibility, 
which is in stark contrast to the experiments that show 
lobe structures for v = 1 as well. 

To summarize, the Keldysh action approach allows us 
to determine the tunneling current through the MZI. The 
theory includes all the effects of the interaction, includ- 
ing charging, renormalization, and dephasing. The de- 
phasing rate grows linearly with voltage; the source of 
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dephasing is the shot noise produced at the first QPC. 
The calculation of this section were carried out under 
the assumption of short-range interaction (which would 
be the case for a system with a nearby metallic gate). Es- 
sential discrepancies between the theoretical findings for 
this model and the experimental data indicate the im- 
portance of long-range (1/r) character of the Coulomb 
interaction. The corresponding model will be studied in 
the next section. 



IV. LONG-RANGE INTERACTION 

In this section, we consider the model of the MZI with 
the long-range Coulomb interaction, 



U (x) 



1 



e Vx 2 + b 2 ' 



(79) 



which is regularized at short distances by the finite width 
of the edge state b. Using the self-consistent electro- 
static picture^ of the quantum Hall edge channels in 
the smooth confining potential, we associate the scale b 
with the typical width of a compressible strip. The latter 
satisfies 



b > a ~ h 



(80) 



where a is the short distance cut-off introduced earlier 
(see also estimates of typical experimental parameters 
in Sec. IIV Cp . The strength of Coulomb interaction is 
quantified by the dimensionless coupling constant 



ve 
e7ri>' 



(81) 



where e is the dielectric constant. As discussed below in 
Sec. IIV CI the edge velocity v (which is also the velocity 
of the neutral mode at filling factor v = 2) is fixed within 
the self-consistent electrostatic picture in such a way that 
a ~ 1. In the momentum space 



2e_ 2 
e 



U (q) = —K (\q\b), 



with the small-5 asymptotic behavior 
2e 2 



Uo(q) 



]n(\q\bd), q<b~ 



(82) 



(83) 



where K is the modified Bessel function of the second 
kind and the numerical constant d — e 7E /2 = 0.89. 

A. Plasmon correlation function 

Consider now the plasmon phase correlation func- 
tion (|2"2"]l . It is given by 




20 40 

t (v/b) 

FIG. 6: Real (upper panel) and imaginary (lower panel) parts 
of the plasmon correlation function plotted versus the di- 
mensionless time at fixed x/b — 50.0; the coupling constant 
a — 1.0 and a/b = 0.1. The dashed line shows the analytical 
approximation 1195 ). 



where we have introduced the phase <f>(q), containing the 
x- and i-dependence, 



(f>{q) = qx-u P (q)t, 
uj P (q) = vq{l + aK (\q\b)). 



(85) 
(86) 



As follows from the integral representation (JHU, & t zero 
temperature Jp(x,t) is an analytic function of time in 
the upper half-plane. We now analyze this function in 
the relevant parameter limit 



i{t,x/v} > b/v > a/v. 



(87) 



The time dependence of real and imaginary parts of 
the phase correlation function Jp(x,t) for typical values 
of parameters is shown in Fig. [51 One observes character- 
istic oscillations which stem from the non-linear plasmon 
dispersion uip(q) = u p (q)q given by Eq. (|86[) . As a re- 
sult, the initially localized density wave packet injected 
into the interferometer arm at some position x' = and 
time t' = is spread after some time t, when it arrives at 
given point x. This is because its constituents — bosonic 
modes, with different wave vectors q — propagate with 
different group velocities. The real part of the phase cor- 
relation function has a pronounced minimum at 
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x/v<t<x/u n (q 2 ) 



FIG. 7: The group velocity of plasmons shown as function of 
the dimensionless momentum for the coupling constant a — 2. 
At qib = 0.59 one has u 9 (qi) = v and for q2b = 1.33 the group 
velocity is at minimum. 



where we have introduced 



d{a) =d e ( Q - 1 )/ Q , 



(89) 



and the numerical constant d has been defined below 
Eq. (j82j) . The time t p (x) is a typical time scale required 
for the fastest plasmon (that has the smallest possible 
momentum q ~ 1/-^) to traverse the distance x. The 
plots at Fig. [S] should be contrasted with the behavior 
of Jp(x, t) in the short-range interaction model given by 
Eq. ([S3")) . In the latter case the real part of Jp(x,t) is a 
smooth function with logarithmic singularity at t = x/v, 
while the imaginary part is piecewise constant. 

The oscillatory behavior of Jp(x,t) can be qualita- 
tively understood using the stationary phase approxima- 
tion to the integral (|84|). For technical details of this 
analysis we refer the reader to Appendix D. The saddle 
point q*(x, t) of the phase <p{q) for given x and t can be 
found numerically as the solution of the equation 



x = u g (q*)t, 
where we introduced the group velocity 
du>p(q) 



dq 

v(l + aK Q (qb) - a (qb)Kx (qb) 



(90) 



(91) 



of the plasmon mode. The plot of the function u g (q) is 
presented in Fig. [JJ As one can see, there are two special 
momenta, q\ and q2, both being independent on cou- 
pling strength a. At q = q\ the group velocity matches 
the drift velocity u g {q) = v, while at q = q 2 the group 
velocity reaches its minimum. Thus, for given x and at 
t < x/v the stationary phase equation (|90[) has a sin- 
gle root (see Fig. [Jj. In the short-time limit one finds 
asymptotically 



q*(x,t) 
(j)(q*) 



bd 



— x/ a. vt 



-aq*vt, t <C x/v. 



(92) 



The unique saddle point leads to approximately single- 
period oscillations in the phase correlation function at 



t p {x) < t < x/v, 



(93) 



as is indeed seen in the numerical plots in Fig. [6] We also 
remark that at t ~ t p (x) the optimal phase <f>(q*) becomes 
of order unity, thus the stationary phase method loses its 
applicability. At still smaller times the integral ([Ml) is 
governed by the contribution of the end point at q = 
and shows no oscillations. On the other hand, in the time 
interval 



x/v < t < x/u g (q 2 ) 



(94) 



one finds (see Fig. [JJ two distinct roots of the stationary 
phase equation (|90[) . In this range the plasmon phase cor- 
relation functions acquires characteristic beatings aris- 
ing from the interference of the two contributions (with 
generally incommensurable phases) corresponding to two 
stationary-point momenta. 

It is shown in Appendix D that at t < u g (q2) the imag- 
inary part of the phase correlation function can be ap- 
proximated in the following way: 



Im Jp~(x,t) = n/2- v^TrRe Erfc^ 



(95) 



where (j> t = 4>{<l*{x,t)) is the stationary phase, Erfc(z) is 
the complementary error function, and the constant 4>o 
is defined as the root of the equation Re [Erfcv / i0o] = 
l/y/2. The above approximation is shown by the dotted 
line in Fig. [6] one finds a good agreement with the exact 
numerical curve. We will use this approximation later 
to evaluate the dephasing action with the logarithmic 
accuracy. 



B. Visibility 

Let us turn now to the analysis of the visibility. For 
that purpose we have to hnd the differential conductance 
in accordance with general relations (foT)]) and (pTTl) of 
Sec. HU We will first perform the analysis in the frame- 
work of the Gaussian approximation (Sec. Ill Cl) . i.e. ne- 
glecting the tunneling action Im A* t in Eq. ([ST]) . Later we 
will include the effects of the non- Gaussian fluctuations. 



1. Filling factor v = 1 

In the case v = 1 the Green function Go at zero bias is 
determined solely by the plasmon part of the total cor- 
relation function J(£), since the free part Jf(0 cancels 
against the bare Green function go'. 



1 



2iri a 



-J>{x,t) 



(96) 
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FIG. 8: Visibility of Aharonov-Bohm oscillations in a sym- 
metric Mach-Zehnder interferometer at v = 1 calculated in 
the model of long-range Coulomb interaction with the cou- 
pling constant a = 1, shown for a/b = 0.1 and different ratios 
of l/b: (l)-30; (2)-50; (3)-100. 



The "lesser" Green function satisfies G'^(x 1 t) = 
(GQ(x,t)) . In the case of finite voltage the above ex- 
pression for the upper arm is modified to 



G%{x,t) = exp<^ -ieV t 



u P (l/L) 



G§(x,t), (97) 



where we took into account the non-equilibrium phase 
shift (j59l- 

In the following we concentrate on the case of symmet- 
ric interferometer, l + = l~ = I. We will show that, in 
contrast to the results of Sec. IIII1 even in this situation 
the visibility is a non-trivial function of voltage. A small 
mismatch in the lengths of upper and lower arms will not 
affect our conclusions essentially. 

In the Gaussian approximation the visibility of inter- 
ference oscillations can be found by integrating Eq. ([6Tj) 
over the time numerically. As discussed above, we neglect 
the dephasing caused by the shot noise at this stage and 
will restore it later. Using Eqs. ([US]) , we obtain for the 
delay time 



At = 



I 



u P (l/L) 



t P (l) . 



(98) 



In a typical experimental layout the distance from source 
to drain L is of the order of MZI size I. Thus, within the 
logarithmic accuracy (i.e. up to a numerical factor in 
the argument of the logarithm) the delay time is actually 
equal to the scale t p {l) defined by Eq. 

The visibility Vg calculated in accordance with the 
definition (|75[) is shown in Fig. [5] for different values of the 
ratio l/b and in Fig. [S] for different values of the coupling 
constant a. In all cases the MZI is symmetric, meaning 
that l + = l~ and both QPCs are equivalent. The plots 
clearly show a non-trivial voltage dependence of visibility, 
with three lobes. The width eo of the central peak scales 
as 



A dip in the visibility at small voltage is the result of a 
zero-bias anomaly in the direct Go(V) and interference 
G$(V) part of conductance. 

As has been already emphasized, the appearance of 
the lobe structure in visibility in the model of long-range 
Coulomb interaction at v = 1 crucially depends on non- 
linearity of plasmon dispersion, which gives rise to the 
oscillatory behavior of the phase correlator Jp(x,t) as 
the function of time. We also note that additional weak 
oscillations in Vg (seen on top of the lobe pattern) stem 
from the beating in the plasmon correlation function in 
the time interval (|94l) . which is due to plasmons with 
short momenta q ~ b . 

Our calculations in this subsection are closely related 
to those of Ref. HH where the effect of non-linearity of 
plasmon dispersion on properties of a v = 1 MZI was an- 
alyzed. The crucial difference is that we also take into ac- 
count the phase shift (|59p representing a non-equilibrium 
charging effect. For this reason, our non-equilibrium 
Green function (p?Tf and the resulting visibility Vg differ 
from their counterparts in the workl28l The voltage de- 
pendence of visibility shown in Figs. [8] and [9] results from 
a combined effect of non-linear dispersion and charging. 



2. Filling factor v — 2 

We have carried out analogous calculations for the 
setup with v = 2. As before, we concentrate on the 
symmetric MZI. In the case of two edge states per inter- 
ferometer arm a tunneled electron excites plasmon and 
neutral modes. The electron propagator at zero voltage 
acquires the form 



G>(x,t) 



1 



-J>(x,t)/2 



2ir\/ iav \Jx/v — t + ia/v 



(100) 




e - 27t/t p (l) 



27rayln(//&) 



(99) 





eV(v/b) 

FIG. 9: Visibility of Aharonov-Bohm oscillations in a sym- 
metric MZI at v = 1 calculated in the model of long-range 
Coulomb interaction for different coupling constant a: (1)- 
2.0; (2)-1.0; (3)-0.5. The other parameters are: l/b = 30 and 
a/b = 0.1. 
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The delay time for the setup of Fig. Q] with one of outer 
channels biased can be found using Eq. (|62l) . 



At 



If 1 1 

2 \u P (l/l) + v 



(101) 



It is instructive to compare this result with Eq. (|98|) : one 
sees a clear manifestation of two modes (plasmon and 
neutral) for v = 2 instead of a single plasmon mode at 
v = 1. The voltage dependence of visibility calculated 
numerically in the Gaussian approximation is shown by 
the upper curve in Fig. 1101 (The lower curve in the same 
plot is the visibility calculated in the instanton approx- 
imation, which takes into account the extra dephasing 
due to the non-equilibrium shot-noise, Sec. IIV B~3l ) One 
can see that the voltage dependence of visibility is qual- 
itatively different from the v = 1 setup. Specifically, the 
visibility Vg(V) shows many oscillations ("lobes") as the 
function of bias with a typical width given by the Thou- 
less energy of the MZI, 



e T h = Try fl 



(102) 



This energy scale differs from the characteristic scale cq 
characterizing lobes in the v = 1 MZI [see Eq. by a 
logarithmic factor. 

The results of this subsection should be also contrasted 
with those of Sec. IIII1 The visibility Vg(V) calculated 
there for the same setup [y = 2) but in the model of 
short-range interaction demonstrates the cos-like oscil- 
lations which, however, do not decay in the case of a 
symmetric MZI. In the 1/r model of Coulomb interac- 
tion considered here, the visibility Vg{V) oscillations do 
decay already in the Gaussian approximation, as seen 
in Fig. 1101 Furthermore, as we are going to show in 
Sec. IIV B 31 the shot-noise dephasing leads in this case 
to an additional suppression of oscillations of visibility 
and to a decay of visibility down to zero at high voltage 
(rather than to a constant as in the Gaussian model). 



3. N on- equilibrium dephasing by the shot noise 

Let us now calculate the suppression of the visibil- 
ity due to non-Gaussian phase fluctuations following the 
method of Sec. HTDl 

In the instanton approximation the optimal phases 
fc ' a are expressed in terms of the phase correlation func- 



tion J(x,t) according to Eq. (|55|). We saw in section HTT1 
that the special properties (|55|) of the causal correlators 
J T / T (x,t) in the case of short-range interaction greatly 
simplify the calculations. The question arises whether 
these properties are also valid in the case of long-range 
interaction. Since J > (x,t) = (J < (x,t))* , we conclude 
that J > and J < may differ only in the imaginary part. 
As is clearly seen in Fig. [6j at x > and t < (or, at 
x < and t > 0) we have ImJp ~ ImJ^ = — ir/2, so 
that the imaginary part of the full correlator, ImJ > , in 
the same range of parameters is almost zero. It follows 




-0,5 0,0 
eV(v/b) 

FIG. 10: Visibility of Aharonov-Bohm oscillations in a sym- 
metric MZI at v = 2 calculated in the model of long-range 
Coulomb interaction with the coupling constant a = 1, shown 
for l/b = 30 and a/b — 0.1. Upper curve: Gaussian approx- 
imation; lower curve: shot-noise dephasing with t1 = 0.2 is 
taken into account. 



that the relations (|65j) are approximately satisfied. This 
is a manifestation of the fact that even in the presence 
of spectral dispersion the edge modes retain their chiral 
character. 

Performing calculation similar to those in Sec. IIII1 we 
arrive again at the result (|70p. which can be cast in the 
form 



Im^* 12 = 2 T 2 1 V{eV)T{l) 1 
where V{eV) is the noise power, 



V(eV) 



di(n<(*) + n>(i) 



and 



T(i) 



r + oo / n 

/ di'sin 2 [ -ImJ> (/,*') -- 
Jo v 



(103) 



(104) 



(105) 



is a characteristic bosonic dwell time related to the size 
of the MZI. In what follows, we will find the dephasing 
action (| 103|) with the logarithmic accuracy (up to a co- 
efficient of order unity under the logarithm). 

First, we evaluate the noise, using its representation in 
terms of interacting Green function 



V{eV) = 4w 2 



dtshi 2 (eVt/2)Re [G>(0,f)] : 



(106) 

To derive this expression, we have used the relation 
G^(0,£) = (G^{Q,t))* and the fact that at zero tem- 
perature and at eV — > the shot noise should vanish. 
To find the interacting Green function, we evaluate the 
momentum integral ((54"]) . At large times t ^> b/v it is 
given by 



J>(t) 



HW 



const . 



(107) 
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Here the momentum q* satisfies the condition | </>(</* )| 
1, which at x = leads to 



(<?*)" 



avt In 



(evt/bdt) 



(108) 



(here e stands for the base of natural logarithms rather 
than for the electron charge). On the other hand, at 
small times t <^ b/v the plasmon correlator is determined 
by the large-momentum behavior of the integrand in Eq. 
(|84p. where interaction is absent and thus Jp{t) should 
coincide with the free result. An approximation for the 
interacting Green function which respects both limits is 
given by 



G>(0,i) 



1 



2ira 1 / v (vt - ia) ln 1/,y (e + ievt/b 



A> (109) 



where we have also taken into account that Gq (0, t) must 
be analytic in the upper half plane of complex t. The 
main contribution to the integral in Eq. (|106l) comes from 
times t ~ 1/eV. Evaluating the integral, we finally find 



V{eV) 



\eV\ 



2Tra 2 /»hi 2/ "(e + v/\eV\bd 



(HO) 



Note that the noise power reveals a weak zero-bias 
anomaly in the chiral ID system, which originates from 
the long-range nature of the Coulomb interaction. 

Let us now evaluate the bosonic dwell time (|105[) . One 
can employ the saddle point approximation (|95[) and 
change the integration variable from time t' to phase <f> 
using the relation (|92|) . We find with the logarithmic 
accuracy 



where C x = 3.96 and C 2 = 2.97. 

Equation ([TUB")) together with Eqs. (fTTU)) and (fTT3)) is 

the main result of this subsection. We see that the non- 
equilibrium dephasing action in the model of long-range 
Coulomb interaction is a sub-linear function of both the 
voltage and the size of the interferometer. Equation (|103p 
tells us that the non-equilibrium dephasing in the MZI 
occurs due to the intrinsic shot noise. This noise is trans- 
ferred by plasmons (in case of v = 1) or by a superpo- 
sition of plasmon and neutral bosonic modes (in case of 
v = 2), which typically propagate through the MZI faster 
than electrons, as is seen from Eq. (|113[) . 

In Fig.[TU]we illustrate the suppression of the visibility 
by the shot noise for the setup with v = 2. The effect of 
shot-noise dephasing leads to suppression of visibility at 
high voltages. However, the lobe structure remains well 
preserved: the dephasing effect becomes strong only for 
voltages much larger than the period of the oscillations 
in visibility. The reason for this is twofold. First, this is 
a factor rf in Eq. (|103[) that was assumed to be small in 
our calculation. (We used the value t? = 0.2 in the plot.) 
We can expect a somewhat stronger dephasing due to the 
shot noise ~ t\ (1 — rf), when the QPC transparency is 
of order 50%. Second, there are logarithmic factors in 
Eqs. (|110|) and (|113[) ; their combined effect at the scale 
eV ~ v/b and at v = 2 leads to the suppression of the 
dephasing action by a factor l/ln 3 (^/6) <C 1. Because 
of this factor, the lobe structure should not be destroyed 
by dephasing (i.e. many oscillations will be seen) even in 
the case of most noisy QPC with transparency of order 
50%, in agreement with experiments. 

A more detailed comparison of our findings with ex- 
periment is presented below. 



I 

a v 



hi 



/ 



b<bd 



(111) 



This enables us to cast the dwell time (|105[) into the form 



T(l) 



I 



x sin 



hr 



/ 



b<j)d 



Re ( Erfcv/i 



(112) 



where <p max ~ a(l/b)/e. This upper limit of integration 
originates from the fact that the integral in Eq. (|105[) is 
actually limited to times t < I /v. Let us now assume that 
the interaction is not too weak, so that </> max 3> 1- This 
is always the case in the experiment, as will be shown 
in Sec. IIV CI Since the integral (|112[) converges at large 
4> (and is dominated by <f> ~ 1), one can now extend 
the integration up to infinity and also neglect the weak 
logarithmic dependence on the phase <j>. In this way we 
find 



(113) 



C. Comparison to experiment 

Let us now discuss the relation of theoretical results of 
Sec. HVI to experimental observations. The striking lobe- 
type structure in the visibility as function of voltage at 
low temperature was discovered for the first time in the 
experiment^. At filling factor v = 1 three lobes were re- 
ported, while at v = 2 five lobes were observed. Later 
experiments of the same group with the use of somewhat 
different layout— showed up to 9 lobes for v = 2. The 
important part of the MZI setup in the experiment^ was 
an additional QPC0, which enabled to apply different 
chemical potentials to the outer and inner edge channels, 
see Fig. [T] These observations were corroborated by sub- 
sequent works^, where an analogous MZI layout with 
QPC0 was studied. 

Our theoretical calculations carried out within the 
model of long-range Coulomb interaction (Sec. IIV[) com- 
pare well with the above experimental findings. On the 
qualitative level, we find three lobes in the case of v = 1 
(Figs. M an d HI) an d structure with many lobes for v = 2 
(Fig. I10[) . in agreement with the experiments. 

Let us also note that a more complicated behavior of 
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the finite-bias visibility was observed in another experi- 
mental workii, where the MZI setup did not contain an 
additional QPCO. The emphasis of this study was on the 
asymmetry of the MZI characteristics when the trans- 
parency of the QPC1 was varied from to 1. We cannot 
treat this feature within our approximation valid in the 
weak-tunneling regime. 

More quantitative comparison of theoretical and ex- 
perimental results requires the knowledge of the edge ve- 
locity v and the range b of the Coulomb interaction. To 
estimate them, we use the results of works 35t39|. In 
Ref. the excitation spectrum of the compressible Hall 
liquid in the classical limit v 1 was found. It consists 
of the magnetoplasmon mode with 

u, (q)=vqln(e-y\q\l), (114) 

and the acoustic spectrum with 



>(g) = vq/n, n > 1, 



where v is the drift velocity, 



v = 



eirh 1 



(115) 



(116) 



and I is the half of the width of the depletion layer be- 
tween the confining gate and the 2DEG. The latter scale 
was found in Ref. |35| and is controlled by the gate voltage, 



(117) 



For the typical concentration no ~ 2 x 10 11 cm -2 , the gate 
potential V g ~1V, and e = 12.5, one thus gets an estimate 
I ~ 110 nm. This result is also believed to be applicable 
in the case of etched mesostructures, e.g. those used 
in the MZI setups. In this case the voltage V g should be 
associated with a work function, which for GaAs-AlGaAs 
heterostructures has the same energy scale ~ 1 eV. 

In the quantum Hall regime, at v = 1 and v = 2, one 
expects that the above modes with n = 0, 1 should match 
the two chiral bosonic modes of our theory. We thus 
see that the drift velocity (|116[) would fix the interaction 
constant a, defined by Eq. ([81]) . to be of order 1. At 
the same time, comparing the dispersion relation (|115[) 
to Eqs. (|82[) and (f56"| . we obtain the estimate for the scale 
b in our theory, b ~ 2e 1 l ~ 400 nm. 

Let us first discuss the case v = 2. For the typical size 
of the MZI, I ~ 10 nm, one gets the ratio l/b ~ 25. This 
parameter controls the overall number of lobes in our 
theory, N ~ I /irb, cf. FigfTUl For sufficiently large trans- 
parency of QPCs the shot-noise dephasing will suppress 
the visibility at sufficiently large voltages corresponding 
to lobe indices > 1/ \n 3 (l/b), see Sec. II V B"3l However, for 
l/b ~ 25 this happens to be not so important for the lobe 
structure. To find the energy scale for the lobe structure, 
we estimate the drift velocity according to Eq. (|116l) that 
yields v ~ 10 5 m/s. This results in the following value 



of the Thouless energy (|102[) : e-rh ~ 20 fieV. This is in- 
deed the typical energy scale seen in the above discussed 
experiments^^. 

We turn to the case of filling factor v = 1. We obtain 
bias dependences of visibility with three lobes in a suffi- 
ciently broad range of parameters (see Figs. [5] and 0, in 
agreement with experimental works^. The energy scale 
for the lobes in our theory is given by Eq. (|99|) , which is 
larger than the Thouless energy exh by a factor ln(Z/6). 
For realistic values of parameters (see above) , this factor 
is ~ 3. In the experiment the energy scales for v = 1 and 
v = 2 are practically equal. This discrepancy may be 
partly explained by the fact that, according to Eq. (|116l) . 
the drift velocity is expected to be twice smaller for v = 1 
compared to v = 2, which reduces the energy scale by 
factor of 2, partly compensating the logarithmic factor. 



V. SUMMARY 

In this paper we have discussed the influence of the 
Coulomb interaction on the quantum coherence in the 
electronic Mach-Zchnder interferometer (MZI) formed by 
integer quantum Hall edge states at filling fractions v = 1 
and v = 2 out of equilibrium. Our main results can be 
summarized as follows. 

1. We have developed the non-equilibrium functional 
bosonization framework which enables us to build 
up the Keldysh action of interacting electrons in the 
MZI. The most non-trivial term in the action is ex- 
pressed in terms of a single-particle time-dependent 
interferometer scattering matrix in the dynamically 
fluctuating field and has a structure of a Fredholm 
determinant similar to those appearing in the the- 
ory of full counting statistics. We have used this 
action to analyze the limit of weak electron tunnel- 
ing between interferometer arms in the case of arbi- 
trarily strong e-e interaction. Our theory contains 
all interaction-induced effects on transport through 
MZI, including charging, dispersion non-linearity, 
and non-equilibrium decoherence. 

2. Restricting at first the theoretical analysis to the 
Gaussian approximation for electron phase fluc- 
tuations, we have readily reproduced the previ- 
ous theoretical results related to oscillations (lobe 
structure) in the dependence of the visibility on 
the bias voltage. Going beyond the Gaussian ap- 
proximation, we have used a real-time instanton 
approach to evaluate the non-equilibrium dephas- 
ing rates that lead to further suppression of the 
Aharanov-Bohm oscillations in conductance with 
the increase of voltage. We have found that the out- 
of-equilibrium dephasing rate is proportional to the 
voltage dependent shot noise of the first quantum 
point contact (QPC1), defining the MZI, and orig- 
inates from the emission of non-equilibrium plas- 
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mons and neutral bosonic modes in course of in- 
elastic electron tunneling. 

3. The results obtained within the model of short- 
range interaction (Sec. IIII[) show strong contradic- 
tions to the experiment. Specifically, in the case of 
equal arms the visibility oscillations in the v = 2 in- 
terferometer do not decay with voltage. For v = 1 
the problem is even more severe, as the visibility 
does not depend on the voltage at all. 

4. Considering the realistic model of the long-range 
(1/r) Coulomb interaction, we are able to explain 
the experimentally observed dependence of the vis- 
ibility of the interference signal at filling fraction 
v = 1 and v = 2. The origin of this effect is found 
to be a combination of three factors: 

(i) the electrostatic phase shift effect, related to 
the charge imbalance on different arms of the 
interferometer, 

(ii) the interaction induced effect of the plasmon 
dispersion, and 

(iii) the out-of-equilibrium decoherence due to the 
intrinsic shot noise. 

Using realistic parameters, we find three lobes in 
the bias dependence of visibility in the case of v = 1 
(Figs. [8] and |9]) and structures with many lobes for 
v = 2 (Fig. I10p . in agreement with the experiments. 
The energy scale for the lobe structure in the visi- 
bility for v = 2 is given by the Thouless energy of 
the MZI which is estimated as ~ 20 /ieV for realistic 
parameters, again in good agreement with experi- 
ment. For v = 1 the energy scale for the lobes 
in our theory is enhanced by a factor ln(Z/6) ~ 3, 
which is partly compensated by a difference in the 
drift velocity (and thus Thouless energy) at v = 1 
and v = 2. There remains some discrepancy in this 
point with the experiment which indicates that the 
energy scales for v = 2 and v = 1 lobe structures 
are practically equal. This issue may be worth fur- 
ther study. 

We anticipate that the approach developed in this work 
will be useful for a much broader class of electronic inter- 
ference setups relevant to current or forthcoming experi- 
ments. The prospects for future research include, in par- 
ticular, interferometers operating in the fractional quan- 
tum Hall regime as well as those built on edge states of 
topological insulators. 
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Appendix A: Regularization of the functional 
determinant 



In this appendix we derive an exact expression for the 
tunneling action At — A— Ao where A is the action given 
by Eq. ([23]) and Aq is its value in the absence of tunnel 
coupling between the interferometer arms. 

Let us denote the first (determinant) term of the action 
(1251 as 



iA c = lndct 



1 + (Ste»*S f - 1)/ 



Then in the absence of tunneling we have 



Ul 0) = lndct 



l + (e 



1)/ 



(Al) 



(A2) 



where 



4>f /h = ~— J o dx'ipf /b (x',t + x'/v F ) (A3) 

are the phases collected along the way from to L with- 
out tunneling. For the full distribution function / one 
can write: 

f{hM) = e^^Mh ~ t 2 y*^\ (A4) 

where <p^(t) = eV ± t, and fp is the Fermi distribution 
function in time domain at zero voltage. Since the gauge 
transformation (|A4j) docs not affect the determinant, we 

can write iA c = IndetM and iAc = IndetMo, with 



M = f > + e lVL S\e l *S s e- lVL f < 
M = f > +e-^ b e i *e^ff < . 



(A5) 



Here /< = fp and /> = 1 — fp, as explicitly defined by 
Eq. (|39[) . Note that after the gauge transformation Mq 
does not depend on voltage. The tunneling action can be 
now represented in the form 



iA t = In det M M n 



(A6) 



and the next step is to invert Mq. 

To perform the inversion, we employ the projection 
properties of the Fermi distribution function at zero tem- 
perature. This approach is well known from the works 
on Fermi edge singularity^ and its recent generalization 
on the matrix case^. In the energy domain, /< is the 
projector on occupied states, whereas /> projects on un- 
occupied states. Therefore, we have /> = />,/< = /<, 
/>/< = = /</>, and /< + /> = 1. The same relations 
hold in the time domain as well, where the product of 
two operators is understood in the sense of convolution. 
Let now ip^ (t) (ifj A (t)) be a function, which is analytic in 
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the lower (upper) complex half-plane. Then the following 
relation hold: 

f<e l ^f< = e^ A /<, />e^ A /> = />e^ A .(A7) 
For instance, let us proof the first relation: 



/<e**7<)(*i,t2) 



1 



dt- 



1 



1 



(2?ri) 2 y h-t + iO t-t 2 + iO 

1 e ^ V ( t2 ) ,v/ , 

f<{t 1 -t 2 )e^ ^. (A8) 



27ri - i 2 + iO 



Due to analytical properties of ^ v (i) this integral is de- 
fined by a single residue at t = t% — iO. 

Furthermore, for any function A(t) with the support 
at real time t, 



A A (t)= / <f<(*-*'M(0 



(A9) 



is the analytic function of complex t in the upper half- 
plane, and 



A^(t) = / dt'Mt-t'W) 



(A10) 

is the analytic function in the lower half plane, respec- 
tively. 

With the help of these zero temperature projection 
properties, we readily obtain the inverse of M . First, 
let us introduce the quantum component of the field tj) 

as 

^9 = $f ~ $b- (All) 

Then we define 

V-' A = /<V^ + x/2, (A12) 
which is analytic in the upper half-plane, and 

i' V = ~f>A - X/2, (A13) 

which is analytic in the lower half-plane. Therefore we 
have 

^ A -V v =^ + X, (A14) 

and thus can write 

M = /> + e-^ V e^ A / < = e-^ 

The inverse of this operator is given by 

Mq 1 = |e-^ V /> + e'^fJ e^ V , (A16) 



e* r f> +e**"f. 



< ■ 
(A15) 



which can be easily checked, using the relations (|A7|) . 
Similarly, we calculate MM^ 1 and find for the tunneling 
action: 



iAt = In dct 



/> 



St&Sft 



' f< 



Let us now introduce the extra gauge phase 

a = f>4>f + f<i>b- 

It enable us to rewrite ip AW in the form 

^ = i>f - A + x/2, 
r = fa -X- x/2. 



(A17) 
(A18) 

(A19) 



Substituting them into Eq. (|A17|) we finally obtain the 
tunneling action in the form (|34[) as stated in the main 
body of the paper. 



Appendix B: Matrix elements of Q 

We define the Q-ma,trix as 

Q = e^Sl( X )S f (x)e-^ = ( * l, V (Bl) 



see Eq. (|4"0|) of the main text. 

To evaluate the action At in the weak tunneling limit 
up to the terms of order of C(r 2 ), one needs to know the 
matrix elements with the following accuracy: 

r = i - tj. - r 2 - Tir 2 jr; r 2 e + r 2 J r{ e j 
+ nr 2 {rfr^e 1 * + r 2 b r{ e - 1 * } e - 1 * (B2) 

+ r 2 2 r 2 h r^e- 4X + rfrfrfe"^ + G(r 3 ), 
T' = in \T\e lx / 2 -r{e~ iX / 2 } (B3) 

+ zT 2 e i$ _ r|e-^ 2 } + 0(r 2 ), 

T = i n {if e-*/ a - rt/e^/ 2 } (B4) 

+ {if e"^ 2 - r^e^ 2 } + C(r 2 ), 

R' = l-r 1 2 -r 2 2 -r 1 r 2 {rJrtV I *+r^r t /e 1 *} 

+ rf I^lfe^ + T 2 T b l r{ f e lx (B5) 

+ nr 2 {r b 2 T\ f e^ + rjr^e"**} e l * + 0(r 3 ). 

They can be easily deduced using the set of rules formu- 
lated in section II. C. 

Appendix C: Tunneling action 

In this Appendix we present technical details of our 
calculations leading to the tunneling action (|43[) . The 
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first step is to evaluate the trace in the expansion (|42|) 
of the functional determinant in the channel and time 
domain using the explicit form of the Q-matrix. A rep- 
resentation of the matrix Q in terms of the "hopping" 
operators Tj defined by Eq. (|4ip is given in Appendix 
B. After straightforward (albeit lengthy) calculations, we 
get 

iA t = J2 r ^ / dt ^ J2 Paa ' ( C1 ) 

ij aa' 

x e- ie ?^n?f(t 1 ,t 2 )e ie ?( t2 \ 

Here P aa ' = 1 if a = a' and P aa ' = — 1 otherwise (a is 
the Keldysh index). Further, 

9t(t) = e+' a (xf,t)-e-' a (xr,t) (C2) 

is the relative phase expressed in terms of the phases 9 ±,a 
accumulated by electron along the paths from scattering 
point xf to the drain at upper and lower arms, 

6 ± ' a {x ± ,t) = -- / dx'<pt(x',t+{x' -x)/v). (C3) 

v Jx± 

The polarization operator n? a has the same structure 
as the operator 11^' , defined by the relations (|44)l , (|45|) 
and pp| . The difference is that fl? a is built up from the 
gauge transformed distribution function f{t, t') 1 given by 
Eq. (|57|) . We also note that the phases 9 ±,a are expressed 
via ip ±,a by the kinematic relation (|C3|) . 

The next step is to show that the action (|C1[) . in fact, 
coincides with its final form (|43[) . For that we use the 
relation 

D* /A (x, t) = ±v- 1 9(±x)6(t - x/v), (C4) 

which follows directly from the definitions (j6]) and the 
standard identity of the Keldysh theory 

D*' A (x,t) = 9(±t) (D>(x,t)-D<(x,t)). (C5) 

The kinematic phase 9 ' a can be now represented in the 
symbolic form as 

9 ± > a (x,t) = (Dfot) (x,t). (C6) 

To proceed further we express the matrix elements of the 
polarization operator ILfj" in terms of the original distri- 
bution functions f^(t,t'), using Eq. This transfor- 
mation reduces the action (|C1[) to the action (|43|) . with 
the phases ±,a equal to 

9 ± ' a (x,i) = \ ± (t-x/v)-9 ± ' a (x,t). (C7) 

Let us now show that the relation (|C7[) is in fact equiva- 
lent to Eq. (0). 

First, using the explicit form (|C4[) of the retarded and 
advanced particle-hole propagator, one can express the 
phases ?pf/ b , defined by Eq. ([36]) . as 

i'b/f = (( D o ~ D o) <pt /f ) (x, t). (C8) 



We also use the fact that at zero temperature the Fermi 
and the Bose distribution functions are closely related to 
each other, namely fp(t— t') = —riB(t — t'). This enables 
us to write the gauge phase (j5B"j) as 

X ± (x-v/t) = ({n B + l)(D£-D«)^){x,t) 

- (n B {D A -D*)yf){x,t). (C9) 

From this relation and Eqs. (|C6[) and (|C7|) one can finally 
see that 

6 ± ^ = -D^f+D<^, (CIO) 

+ (Gil) 

which agrees with the desired relation To obtain 

these identities, we have used the standard properties of 
the Keldysh propagators, 

D< = n B (D§ - D£), D> = (n B + - D$), 

Dl = D* + D<, = Dq — Dq. (C12) 

This completes our proof of equivalence between two 
forms of the action, Eqs. (|C1|) and (|43l) . 

Appendix D: Stationary phase method 

To describe the oscillations of the plasmon correlation 
function Jp(x, t), we analyze the integral ([54")) in the 
short time limit t <C x/v. Then the optimal momentum 
is small, <C b , which enables us to write 

<f)(q) = 4>*g(q/q*), (Dl) 

where ^> 1 is the stationary phase and 

fl (A) = A(l-lnA). (D2) 

Let us consider 

d,j, t J> = -i dX(l-\n\)e^ g ^e- aq * x . (D3) 
Jo 

The main contributions to this oscillatory integral come 
from the stationary points and from the end points. 
Here we concentrate on the stationary phase contribu- 
tion, which is responsible for the oscillations in the cor- 
relation function. Thus we get 



fy.J> = - t y£^-e*.. (D4) 

Integrating back this equation and taking into account 
the boundary condition, Im Jp — > ir/2 for t — > oo (0* — > 
— oo), one arrives at 

Im j£ = | - V2ttRc Erfc( , (D5) 
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where Erfc(z) is the complementary error function. This 
formula does not apply at t — > (0* — > 0), where the 
stationary phase is small and the associated approxima- 
tion breaks down. As is seen in Fig. [BJ at <^>* ~ 1 the 
oscillatory behavior crosses over into a smooth one, and 
with further decrease of t the function ImJp saturates 
at — 7r/2. Wc can approximate this crossover at 0* ~ 1 
and saturation at </>* ~ by replacing <f>* in Eq. (|D5[) by 
0* — 0o i where the phase 0o ~ 0-12 satisfies the equation 

ReErfc(v / #o) = 1/V2. (D6) 



To improve further the accuracy of the approximation, 
one can also replace the limiting expression for the sta- 
tionary phase (|92|) (which was found using the logarith- 
mic approximation to the plasmon dispersion relation) 
by its exact value found from the numerical solution of 
Eq. (|90]l. This yields the final approximation, which is 
shown in Fig|6]by the dotted line. 
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